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ABSTRACT 

Faraday rotation is a valuable tool for detecting magnetic fields. Here the technique is con- 
sidered in relation to wind-blow bubbles. In the context of spherical winds with azimuthal or 
split monopole stellar magnetic field geometries, we derive maps of the distribution of position 
angle (PA) rotation of linearly polarized radiation across projected bubbles. We show that the 
morphology of maps for split monopole fields are distinct from those produced by the toroidal 
field topology; however, the toroidal case is the one most likely to be detectable because of its 
slower decline in field strength with distance from the star. We also consider the important case 
of a bubble with a spherical sub-volume that is field-free to approximate crudely a "swept-up" 
wind interaction between a fast wind (or possibly a supernova ejecta shell) overtaking a slower 
magnetized wind from a prior state of stellar evolution. With an azimuthal field, the resultant 
PA map displays two arc-like features of opposite rotation measure, similar to observations of 
the supernova remnant G296. 5+10.0. We illustrate how PA maps can be used to disentangle 
Faraday rotation contributions made by the interstellar medium versus the bubble. Although 
our models involve simplifying assumptions, their consideration leads to a number of general 
robust conclusions for use in the analysis of radio mapping datasets. 

Subject headings: ISM: supernovae remnants - Polarization - Radio Continuum: stars Stars: circum- 
stellar matter - Stars: magnetic fields - Stars: winds, outflows - 



1. Introduction 

Magnetism plays an important role in the lives 
of stars, frequently in the form of its influence on 
angular momentum transport for star formation 
and evolution (e.g., Ghosh & Lamb 1979; Bland- 
ford & Payne 1982; Bodenheimer 1995; McKee & 
Ostriker 2007; Meynet, Eggenberger, & Maeder 
2011; Matt et al. 2012), and also in terms of 
hot plasma generation (e.g., Davidson & Ostriker 
1973; Babel & Montemerle 1997ab; Townsend, 
Owocki, & ud-Doula 2007; Li et al. 2008; Giidel 
& Naze 2009). Of interest to this paper is the 
growing body of evidence for magnetism among 



massive stars. Direct detections of magnetism in 
normal stars (i.e., not compact objects) are in 
large part relegated to measuring circular polar- 
izations in spectral lines arising from the Zeeman 
effect (see Donati & Landstreet 2009). The first 
detection of magnetism in a star besides the Sun 
dates back to Babcock (1947). Since then, the 
field has exploded. A recent review of the current 
state of the subject for non-degenerate stars can be 
found in Donati & Landstreet (2009). The key re- 
sult now is that whereas magnetic detections were 
mainly limited to stars with surface fields in the 
kilo-Gauss range, modern instrumentation and di- 
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agnostic strategies regularly achieve highly signif- 
icant detections in the regime of 100's of Gauss 
(e.g., Donati & Collier Cameron 1997), and in 
some special cases much less (e.g., Scnnhauser & 
Berdyugina 2011). 

Numerous direct detections of surface mag- 
netism in massive stars have been reported, with 
some recent examples being Alecian et al. (2011), 
Petit et al. (2011), Scholler et al. (2011), Hubrig 
et al. (2011, 2012) Grunhut et al. (2012a, 2012b, 
2012c), Wade et al. (2011; 2012). These suc- 
cesses have correspondingly motivated theoretical 
studies to understand the origin of these fields 
for massive stars (e.g., MacGregor & Cassinclli 
2003; Braithwaite 2006; Cantiello et al. 2009), 
their influence on massive star evolution (Maeder 
& Meynet 2003; Yoon, Dicrks, & Langer 2012), 
and connection to other observational phenom- 
ena such as X-ray emissions (Babel & Montmerle 
1997b; Gagne et al. 1997; Gagne et al. 2005; Fa- 
vata et al. 2009; Ignace et al. 2010; Oskinova et al. 
2011; Gagne et al. 2011; Wade et al. 2012; Grun- 
hut et al. 2012; Ignace, Oskinova, & Massa 2012) 
and aspherical wind flow (Poe, Friend, & Cassinclli 
1989; Shore & Brown 1990; Chevalier & Luo 1994; 
ud-Doula & Owocki 2002; Townsend & Owocki 
2005; Brown, Cassinelli, & Maheswaran 2008; ud- 
Doula, Owocki, & Townsend 2008; ud-Doula et al. 
2012). 

There are other diagnostics of stellar mag- 
netism that have or could complement the Zeeman- 
based approach. Non-thermal radio emissions 
from massive star colliding wind binaries have 
been used to infer stellar magnetism (e.g., Williams 
et al. 1997; Dougherty & Williams 2000; De 
Becker et al. 2006; van Loo, Runacres, & Blommc 
2006). The Hanle effect is a weak Zeeman effect 
pertaining to the influence of magnetic fields on 
the linear polarization in spectral lines. The ef- 
fect is sensitive to magnetic fields in the 1-100 G 
range (depending on the Einstein A-value of the 
lines being measured) and has been successfully 
used in studies of solar magnetism for decades 
(e.g., Sahal-Brcchot, Bommier, & Leroy 1977; 
Stenflo 1982; Berdyugina & Fluri 2004; Trujillo 
Bueno et al. 2005). There is a small but grow- 
ing literature on its potential application to other 
stars (Ignace, Nordsieck, & Cassinelli 1997; Lopez 
Ariste, Asensio Ramos, Gonzalez Fernandez 2011; 
Ignace et al. 2011; Bommier 2012; Manso Sainz & 



Martinez Gonzalez 2012). 

Another important method for measuring as- 
trophysical magnetic fields is Faraday rotation. 
The effect refers to how the line-of-sight (LOS) 
magnetic field component rotates the position an- 
gle of linear polarization for a beam of radiation. 
The amount of rotation is also proportional both 
to the electron density (hence operates only in a 
plasma) and to the path length through such re- 
gions. Importantly, the amount of position an- 
gle (PA) rotation scales with the square of the 
wavelength of observation, A 2 . In this way most 
applications measure the polarization PA for a 
range of wavelengths to derive a quantity called 
the "rotation measure" (or RM; see the following 
section) that encodes information about the inte- 
grated product of the LOS field component and 
electron number density. 

Most applications of Faraday rotation are for 
interstellar or extragalactic studies. Space pre- 
cludes a comprehensive review of this literature; 
discussion of the state of the field, with refer- 
ences therein, can be found for example in Car- 
illi & Taylor (2002), Han et al. (2006), and Beck 
(2012). Here attention is focused on the potential 
of Faraday rotation as a probe of stellar magnetism 
in wind-blown bubbles. There has been several 
recent observational developments that speak to 
the interaction of stellar winds or supernova ex- 
plosions with the surrounding interstellar environ- 
ment and the capacity of probing these interac- 
tions with Faraday rotation. 

Ransom et al. (2008, 2010) have conducted 
studies of Faraday rotation effects arising from 
planetary nebulae (PNe). The relative motion 
between a PN and the surrounding interstellar 
medium (ISM) alters the strength and direction 
of the interstellar magnetic field, leading to vari- 
ations of the polarization PA across the PNe and 
the trailing tail it creates in relation to the PAs 
surrounding the structure. Although implicitly 
placing constraints on the stellar magnetic field 
interior to the nebula, the PNe act essentially as 
perturbers, on the scale of a few parsecs, of the 
local medium, with Faraday rotation serving as a 
probe of the resultant disturbances. 

Regarding massive star influences, Savage, 
Spangler, & Fischer (2012) report on an exten- 
sive study of Faraday rotation for the Hn region, 
the Rosette Nebula. In this application a cluster 
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of massive stars lead to a wind-blown bubble and 
surrounding photoionized region by the central 
OB stars. In contrast to the PN studies of Ran- 
som et al. (2008, 2010), whose analysis was based 
on PA maps with the diffuse Galactic synchrotron 
used as a source of linearly polarized radiation, the 
Rosette study of Savage et al. analyzed data from 
an array of sightlines to background extragalactic 
sources of polarized emission that intercept the 
nebula and surrounding region. Their data are 
consistent with PA changes across the nebula, on 
a scale of ~ 10 pes, arising from the presence 
of the bubble and its impact on the interstellar 
magnetic field. 

For the case of Faraday rotation as a probe 
of stellar magnetism, Harvey-Smith et al. (2010) 
discuss an antisymmetric RM morphology across 
the supernova remnant (SNR) G296. 5+10.0. The 
SNR nebula has two prominent emission arcs on 
opposite sides of a symmetry axis, yet the PA ro- 
tations of either arc are oppositely oriented. This 
kind of pattern would be expected from a mag- 
netic field that reverses its LOS polarity from one 
arc to the other. The authors associated the pat- 
tern with the magnetic field of a slow magnetized 
red supergiant wind. The ejecta shell from the SN 
has "swept up" the supergiant wind into the ob- 
served shell. To reproduce the RM pattern, the 
stellar magnetic field would have to be toroidal to 
produce the observed polarity change between the 
emission arcs. The authors derive an expression to 
relate the observed scale of PA rotations to proper- 
ties of the shell and the stellar wind. The observed 
change in amplitude of the RM by approximately 
40 rad/m 2 indicates that a surface stellar magnetic 
field on the order of 500 G could account for the 
observations. 

In this paper we explore further the idea of 
Faraday rotation as a means of deriving infor- 
mation about stellar magnetism in the large-scale 
wind that might sometimes be observed in wind- 
blown bubbles and SNRs. In section [2] a brief 
review of the expressions describing Faraday ro- 
tation are given. A derivation of PA rotation 
maps for an ionized and spherical stellar wind 
with azimuthal (or toroidal) magnetic fields is pre- 
sented in section [2~T1 PA rotation maps for a split 
monopole are presented in section 12.21 as a con- 
trast case. Returning to the azimuthal field in sec- 
tion I2.3[ simple insertion of a spherical subvolumc 



taken to have no magnetic field is used to simulate 
a two-wind interaction like that of a SNR or any 
scenario in which a fast flow overtakes a slower 
one from an earlier stage of stellar evolution. In 
section [3] applications of our results are discussed, 
with concluding remarks given in section [4] 

2. Theoretical Models 

The standard expression to represent the po- 
larization PA rotation, ip, arising from Faraday 
rotation along a sightline is 

i> = Tpo + RM X (-^\ , (1) 

where ipo is the orientation of a background po- 
larization upon which Faraday rotation operates, 
and RM is the "rotation measure" given by the 
path length integral of 

EM.O.Slradm-/ (*|) (-^) , 

(2) 

where By is the LOS component of the magnetic 
field, n c is the electron density, and path length. 

The rationale for this particular formulation is 
that observers are frequently interested in deriv- 
ing RM from the radio data as a constraint on the 
magnetic field strength (modulo the electron den- 
sity, whose value may be constrained from other 
considerations). One measures values of ifi for a 
range of wavelengths; plotting tp against A 2 in a 
log-log plot should then yield a straight line whose 
slope is the rotation measure, RM. If RM > 0, 
then the net Faraday rotation from all of the par- 
allel field components along the LOS, some be- 
ing positive and some being negative, in regions 
where there are electrons has resulted in a counter- 
clockwise rotation of the background linearly po- 
larized radiation. If RM < 0, the net effect is a 
clockwise rotation. The sign of ip, or alternatively 
RM, is sensitive to the electron-density averaged 
net LOS field component along a given sightline. 

The distinction being drawn here is that our 
models emphasize differential Faraday rotation 
(e.g., Minter & Spangler 1996). This amounts to a 
map of how RM varies between neighboring sight- 
lines, which is equivalent to analyzing maps of ip 
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at a fixed wavelength. Consequently, it is conve- 
nient for our purposes to reformulate the effect of 
Faraday rotation as: 

*- 555 /(!)(£)* < 3 » 

where no and Bq are chosen normalizations, zq(X) 
is a relevant length scale for the problem, dz is in 
the same units of zq, and the factor of it is made 
explicit here both to indicate that tp is measured in 
radians and to highlight the fact that polarization 
PA is degenerate in multiples of it (or 180°). 

It is straightforward to convert such a map in 
PA to one in RM. If one evaluates i/j(\o) at a ref- 
erence wavelength Ao, then 

yj(X) = V(A ) x * (4) 

and so 



The following discussion of Faraday rotation in 
stellar bubbles adopts the use of equation ([3]). 

We focus on three illustrative cases to high- 
light expected PA map morphologies as motivated 
by current observations. First the case of an az- 
imuthal field distribution is considered. Since the 
field strength diminishes only as r _1 , it is this 
component that is most likely to have observa- 
tional relevance at the 1-10 parsec scales of stellar 
bubbles. Second, as a contrast case, PA maps are 
derived for a split monopole field. This case leads 
to results that are morphologically distinct from 
the azimuthal field case. The field drops much 
faster with distance from the star, as r~ 2 , mak- 
ing this case essentially unobservable; still, the 
split monopole has value in providing insight into 
the range of morphological possibilities. Finally, 
we consider again the azimuthal field scenario but 
now imposing a central "cavity" that does not con- 
tribute to Faraday rotation (e.g., a central region 
devoid of a magnetic field). This approximates 
a "swept up" field from a wind-wind interaction, 
such as may occur for a SNR or a PN, where only 
the wind from an earlier phase of stellar evolution 
has a relevant magnetic field. 

We stress that spherical symmetry is adopted 
for the bubble shape and density distribution. It 



is not the goal of this paper to reproduce the ob- 
servations for any particular object. A spherically 
symmetric density profile provides a "controlled" 
environment for which to evaluate and gain in- 
sight into Faraday rotation through stellar bub- 
bles. Real situations may involve a broad range 
of additional (and potentially important) consid- 
erations, such as: aspherical density distributions 
(e.g., bipolar flows or clumping); radius-dependent 
and/or aspherical ionization effects; correlated be- 
havior between density, ionization, and/or a dy- 
namically relevant magnetic field (e.g., how a 
toroidal field can lead to axisymmetric bubbles as 
described in Chevalier & Luo 1994). Our spherical 
results provide a suite of baseline cases that are an- 
alytic or semi-analytic in which one can evaluate 
modifications to predicted PA maps arising from 
these more complex factors. 

Before proceeding it is useful first to review the 
different contributing components to Faraday ro- 
tation that affect the final observed PA at the 
Earth. The underlying assumption is that Faraday 
rotation acts to rotate the PA of linearly polarized 
radiation as it passes through magnetized and ion- 
ized regions. The background source (whether a 
point source or diffuse synchrotron background) 
has some initial position angle V'orig- This value 
receives an additional contribution ipiSM owing to 
the ISM. The total PA measured at the Earth is 
then 

Iproeaa = V'orig + IpISM, (6) 

where both f/Vig and ipiSM are independent quan- 
tities that can vary from one sightline to the next. 

Now for sightlines that intersect a stellar bub- 
ble, one must subtract the ISM contribution for 
that segment of the sightline passing through the 
bubble, which we signify as StpisM- Then the con- 
tribution made by the bubble, V'bub: must then be 
added. The total PA at the Earth now becomes: 



V'meas = V'orig + feilf ~ ^IpISM + V-'bub- (7) 

The calculations in the following sections present 
maps specifically of V'bub- Section [3] deals with the 
fact that the bubble contribution alone is not what 
is actually measured. 
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2.1. Azimuthal Fields 

Modeling the PA morphologies of stellar bub- 
bles requires the introduction of two coordinate 
systems: the observer system and the stellar one. 
For the observer Cartesian, cylindrical, and spher- 
ical coordinates are adopted as (x, y, z), (w, a, z), 
and (r, 6, a), where the origin is the bubble center. 
Here z is the observer axis, with the Earth at +00. 
The angle a is measured counter-clockwise about 
this z-axis. Then 9 is the polar angle from the ob- 
server axis. As we will be discussing sightlines in- 
tersecting the stellar bubble, the cylindrical radius 
m will be the impact parameter for such rays. For 
the star the corresponding Cartesian and spherical 
coordinates are (x*, j/,, z*) and (r, tp). 

For the transformation between these coordi- 
nates, we choose y = y*. Using unit vectors, the 
viewing inclination angle i between the z and z* 
axes is given by 

cosi = z-z„. (8) 

Transformations between the angular quantities 
can be obtained with spherical trigonometric re- 
lations, that will be used as needed. 

In this first example, an azimuthal stellar mag- 
netic field is considered. The vector field is B = 
B v tp. Calculation of Faraday rotation along a 
sightline requires determination of the LOS field 
component, which is given by 



& Cassinelli (1998) based on wind compression 
theory (Bjorkman & Cassinelli 1993). The Ig- 
nace et al. model assumes a magnetic field that 
is dynamically negligible as compared to the wind 
flow. Assuming flux freezing, the axisymmetric 
field topology can be derived, a result they refer 
to as "WCFields". Their model includes a pa- 
rameter for the axisymmetric density distribution. 
For our purposes this parameter, d/i/d^o, is set to 
unity for a spherical wind; then using their asymp- 
totic formula (eq. [21] of Ignace et al. 1998), the 
azimuthal field becomes simply 

B v (r,0) =5eff^ sintf, (12) 

where 1? is the co-latitude on the star, signifying 
that the maximum toroidal field strength occurs 
at the rotational equator of the star; and for the 
model of Ignace et al. , the conveniently defined ef- 
fective surface field strength, B e g, depends on the 
actual surface field strength the stellar rota- 
tion speed w rot , and the wind terminal speed Wqc, 
with 

B cS = B, . (13) 

Calculation of the Faraday rotation proceeds 
from a sightline-depcndcnt integration through 
the bubble, as given by 



B ■ z = B v (tp ■ z). 



(9) 



The transformation between the observer and stel- 
lar coordinates is needed to evaluate the preceding 
dot product; the rotation matrix between coordi- 
nate systems is given by 



^bub(ro,aO 



7T 



n c (r) 



(z-tp) dz, 
(14) 

with n e = n w (i?»/r) 2 , where n w is a density scale 
associated with the wind. This integral reduces to 



y* 
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sinz \ 







1 






— sinz 





cos I I 





Consequently, one obtains 



S|| = — B v sin tp sin i. 



(10) 



(11) 



In addition to the LOS field component, the 
distribution of the field strength with location 
about the star is also needed. We adopt the kine- 
matic prescription put forth in Ignace, Bjorkman, 



V'bub(n7,a) = — — f (— ) sin 1? sin tp sin idz. 
Zb J \ r J 

(15) 

Note that sin$sin<p = y*/r = y/r for our coor- 
dinate system definitions. Given that y _L z, the 
coordinate y can be factored out of the integral 
which now becomes 



V>bub(ro,a) = -7r sini — 



Zb 
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To solve this equation, we note that r 2 = vo 2 + z 2 , 
where w is the impact parameter of the sightline 
under consideration, and therefore a constant of 
the integration. The sightline enters the bubble at 
+zq and exits at — Zq, where zq(vj) = \J R 2 — vj 2 
(see Fig. [IJ. Back- front symmetry of the integra- 
tion gives 



V'bubCn?, a) = — 27r sini — / 
z b Jo 



2o(ro) / R 2 V dz 



m 2 + z 2 J R*' 
(17) 

which, after some rearrangement, has the solution 



ipbnb(x,y) 



-2tt 



R* 



(ir/4-9 /2 



1 



R* 
w 

sin 2^o 



(18) 



where tan#o = I z§(vj) and w = y 1 x 2 + y 2 . 

There are several key comments to be made 
about this solution. 

• First, the appearance of the factor of y 
means that the PA map is left-right anti- 
symmetric about the line of x — in the 
plane of the sky. 

• Second, the overall morphology of the PA 
map is independent of viewing inclination. 
The inclination angle appears in the solution 
only in the multiplicative factor sini, acting 
as an amplitude scale. Consequently, the PA 
map for the edge-on view is the same map 
that results for any other inclination, just 
the amount of PA rotation is reduced for ev- 
ery sightline by sini. 

Figure [5] shows a false-color plot of the PA map 
(or after proper normalization, the RM map) for 
the solution of equation (|T%1) . Based on our con- 
ventions, the positive :r-axis points down in this 
figure, and the positive y-axis points right; to- 
gether these give positive z toward the observer. 
The map is indeed left-right antisymmetric and 
top-down symmetric. Additionally, because the 
azimuthal field is everywhere perpendicular to the 
LOS for x = 0, V'bub = along the vertical that 
passes through the map center. 



No absolute scale is given in this figure as it 
depends on a number of wind and star parameters, 
such as the surface field strength, the mass-loss 
rate, wind speed, and stellar radius as well as the 
bubble radius (hence, its age). We return to the 
expected level of PA and applications in section [3] 

2.2. Split Monopole Fields 

It is useful to consider a different field topology 
to explore the range of PA map morphologies that 
can result in stellar bubbles. An azimuthal field 
exhibits the most shallow radial decline expected 
for a field that is carried out by a stellar wind. 
The next most shallow decline would be a split 
monopole. This is a radial magnetic field with a 
field strength B x = B* (i?*/r) 2 , but the polarity 
changes sign from one hemisphere to the other. So 
the field is outward directed in one hemisphere, 
but inward in the other. Such a magnetic con- 
figuration is what would be expected of a strong 
stellar wind that distorts a dipole field at the star 
into a radial geometry (e.g., ud-Doula & Owocki 
2002). 

Calculation of PA maps proceeds as before, ex- 
cept that now there are different factors appearing 
in the integrand because of the new field topology. 
Instead of z ■ (p for the azimuthal field case, we now 
have z-f = z/r for a radial field. The integral takes 
on the form: 



bub 



Zhub 



n c (r) 



Br(r) 



Z ~) dz. 

r / 



(19) 

Before proceeding, the split monopole case of- 
fers a new wrinkle for the calculation of the Fara- 
day rotation. The radial field switches polarity 
between hemispheres. Consequently, at a general 
viewing inclination, the net Faraday rotation can- 
cels identically for sightlines which do not inter- 
cept the magnetic equator. Such sightlines en- 
ter and exit the bubble in a hemispherical cap of 
just one field polarity. As a result, there is as 
much clockwise PA rotation through, say, the first 
half of the pathlength as there is counterclockwise 
contribution through the second half. The radial 
field ensures that for every value of B« along the 
path, there is a corresponding value of —B» at a 
reflected position in a back-front sense along the 
path. These two positions of opposed LOS field 
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components occur at the same radius from the 
star, and therefore occur at the same density. A 
polarity switch, and therefore a net Faraday rota- 
tion, only occurs for sightlines that intercept the 
magnetic equator. 

The result of all this is that a pair of truncated 
hemispherical caps appear at top and bottom in 
the PA maps for the split monopole case. The ex- 
tent of these truncations depends on the viewing 
inclination. For the pole-on case, every sightlinc 
passes through the magnetic equator, and there 
are no truncation zones. The extreme opposite 
case is the edge-on view; here, no sightlines pass 
through the magnetic equator, and so i/'bub is iden- 
tically zero everywhere. In observer coordinates 
the truncation occurs for x > ±Rcosi. 

For sightlines that do intercept the magnetic 
equator, the location along the path where this 
occurs, z oq , must be known so that the sign change 
is properly taken account in the integration. For 
sightlines with \x\ < -Rcosi, the integral becomes 



"0bub = -R* 



r+z (ro) 



— dz 



zo(ro) 




where a sightline with impact parameter w at ob- 
server azimuth a for a split monopole viewed at 
inclination i will intercept the equatorial plane at 
location 



z cq = —vo tanz cos a. 



(21) 

Evaluating the integral of equation (I2TJ1) . along 
with some algebraic manipulation, leads to 



^ b b - -— (—] (—) \( R s 

3 \ z b / \ R/ \x 2 +y 2 cos 2 i 

(22) 

This solution is displayed in Figure [3] at a viewing 
inclination angle of i — 60 c irc. 

Unlike the case of an azimuthal field, the scal- 
ing of the results for the PA is more complicated 
with viewing inclination. The truncated caps are 
inclination dependent. The PA maps for a split 
monopole are markedly different from that of an 
azimuthal field as seen in Figure [5J Whereas 
the azimuthal field produces an antisymmetric 
morphology, all the Faraday rotations for a split 



monopole are of the same sense: either everywhere 
clockwise or everywhere counterclockwise. 

As noted previously, the effects of Faraday rota- 
tion are much smaller across the bubble for a split 
monopole than for an azimuthal field. However, 
assuming that ^bub ~ for the vast majority of 
sightlines in the split monopole case, the bubble 
still influences the PA map of the region because 
of differential RM effects for sightlines intercept- 
ing the bubble versus those that do not. In other 
words even if ipbub <S 1, the mere presence of the 
bubble means that Si^ism may still be significant. 
More will be said about this in section [3] 

2.3. "Swept Up" Fields 

For the last application, consider again the az- 
imuthal field case. To simulate a two- wind inter- 
action, where an inner wind (or a supernova ejecta 
shell) is blown into, and sweeps up, an outer wind 
from a earlier stage of stellar evolution, we impose 
a central spherical region that does not contribute 
to Faraday rotation. This could be either because 
the central region has no magnetic field or because 
it consists of neutral gas (the latter being unlikely 
for applications to massive stars). The overall sce- 
nario is intended to approximate situations like the 
production of a PN, or a fast wind from a blue su- 
pergiant overtaking a slower denser red giant wind 
(e.g., Chita et at 2008), or a SN explosion that ex- 
pands into the wind of a progenitor phase. 

To be specific, we define an interior boundary 
of radius R\ < i?2, now with i?2 the outer bub- 
ble radius. For r < R\ the region is effectively a 
"cavity" in terms of its contribution to the Fara- 
day rotation along a sightline. The polarization 
PA, V'bub) n ow becomes 

^bub oc y sin i x / dz/r 3 , (23) 

J zi 

where z\ = for vj > R\ 1 but for zu < R\, z\ 
takes on the value 



R\ 



(24) 



Detailed steps for the derivation of V'bub as a func- 
tion of sightline, along with a generalization to a 
power-law density, are given in the Appendix. Fig- 
ure 2] shows examples of the resulting PA maps for 
cavities of different relative extents with Ri/R 2 = 
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0.25, 0.5, and 0.75. The maps remain left-right an- 
tisymmetric as in section [2~Tl but now two maxima 
appear in the form of arcs. 

Not surprisingly, these maps are quite reminis- 
cent of the morphologies seen in SNR G296. 5+10.0 
(Harvey-Smith et al. 2010). In that report a 
discussion and derivation similar to this paper 
are presented. Harvey-Smith et al. introduced a 
toroidal held to explain the antisymmetric depen- 
dence of RM observed in the SNR. Those authors 
additionally included density and magnetic held 
enhancements arising from a shock. These are ad- 
ditional relevant ingredients that we have not in- 
corporated here, but such effects could be included 
either in a phenomenological way or using detailed 
dynamical simulations. Our approach highlights 
the robust nature of the antisymmetric PA map 
morphology provided by an azimuthal held along 
with the appearance of arc-like features arising 
from the presence of a "cavity" , a region devoid 
of any contribution to Faraday rotation. 

3. Discussion 

The objective of this paper has been to elu- 
cidate the effects of Faraday rotation through a 
stellar bubble by evaluating of PA maps (or equiv- 
alently RM maps) under certain simplifying as- 
sumptions to focus on broad morphological trends. 
As such, spherical symmetry has been assumed for 
the geometry of the bubble and the run of density, 
specifically a wind type density that drops as r~ 2 
with radius. At this point it is worth commenting 
on the limitations of this approach, the kinds of 
effects that need to be included for application to 
particular types of bubbles, and strategies for best 
extracting information about the properties of the 
bubble and its immediate environment. 

The case of an azimuthal stellar magnetic held 
is likely the only one of observational relevance be- 
cause of its more gradual decline in field strength 
as r -1 to allow for detection at the 1-10 pc scale 
of stellar bubbles. A robust prediction of even 
the simplified models presented here is that a sign 
reversal in polarization PA (or RM) is expected 
when an azimuthal field plays a role in the Fara- 
day rotation. However, that sign reversal will only 
be achieved if the f/'bub term can be isolated. 

Imagine a situation in which the ISM along a 
sightline passing near to a bubble, but without 



intercepting the bubble, produces a PA value of 
the form: 

IpISM =tt x L /zism, (25) 

where Lq is the relevant length scale over which 
one can define an average product of the interstel- 
lar magnetic field and electron number density. In 
other words: 

(B ]{ n e ) = -3- / B\\n c dz. (26) 
L o Jo 

Assuming that this average value is what would 
have been sampled for a sightline through the bub- 
ble, and further assuming spherical symmetry, the 
decrement to the interstellar PA rotation caused 
by the presence of the bubble as a function of im- 
pact parameter is 



Sil>isM(™)=fl>iSM-j-yl-^. (27) 

The ISM is highly inhomogeneous in terms of 
its magnetic field distribution, density, and ioniza- 
tion. Consequently, one hardly expects that n c B^ 
at the location of the wind bubble would equal the 
pathlength-averaged value (B\\ n c ). Still, it is use- 
ful to consider the resultant PA map from the com- 
bination of the different contributions by adopting 
such a scenario. Choosing -0-TSAilll = 1000° and 
2R/Lq = 0.01, Figure [S£i shows a PA map and 
contours for an azimuthal field case, similar to the 
style of previous figures, but now with the Faraday 
rotation contribution from the ISM included. Ad- 
ditionally, Figure \Sjp displays what the measured 
PA would be along an axis at x = 0.1R. 

Figure [5] is important for illustrating how to in- 
fer the various components to the measured value 
of V'mcas- If one could measure the PA map across 
a wind-blown bubble with an azimuthal field, the 
following outlines how to decompose the trace of 
the signal in Figure[5}) to determine ipisM, StpisM, 
and ipbub- Note that it is assumed that the inter- 
stellar field and electron density in the vicinity of 
the bubble is approximately constant (or at least 
smoothly varying). 



1 Recall that polarization PA is degenerate in multiples of 
180°; here tp IS M = 1000° amounts to a PA of 100°, as 
seen in Fig. [5p. 
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• For \y\ > R, the PA is set by sightlines 
that do not intercept the bubble and there- 
fore sample the interstellar field and electron 
density. 

• For \y\ < R, the sightlines intercept the bub- 
ble. The azimuthal field leads to ^bub that 
is predicted to be antisymmetric, but the 
decrement term StpisM is symmetric. The 
result is a left-right asymmetric profile. 

• The two peaks in Figure [5jb) can be used to 
deduce the amplitude of V'bub- The ver- 
tical blue double-arrowed line represents 
the total change in the PA rotation ow- 
ing to the bubble contribution. Introduc- 
ing ip max (x) = |^bub(a;, ?/max)| as being the 
maximum PA rotation through the bubble 
at location (x, ±y max ), the blue arrow has 
an amplitude of 2?/> max . The value of V'max is 
related to the scale length z\, that depends 
on the wavelength of observation, the den- 
sity scale of the bubble, and importantly the 
surface field strength of the star. Bisection 
of the vertical arrow, as illustrated with the 
horizontal dotted line, yields the value of 

• The bisection mentioned in the previous 
point gives the maximum decrement due 
to 5-ipisM- With spherical symmetry, one 
should expect a distribution of this decre- 
ment as indicated by the green dashed line. 
Bear in mind that the particular example 
of Figure [5] implicitly assumes a local inter- 
stellar field of positive polarity, which leads 
to the bowl-shaped decrement as shown. If 
the local field were of the opposite polarity, 
the "decrement" would in fact produce an 
inverted bowl shape, meaning that the green 
line for SipisM would lie above ipiSM- Con- 
sequently, the decomposition process gives 
both the strength of the local magnetic field 
and its polarity. As one last comment, it is 
possible that no decrement is found. Such a 
result would arise if the interstellar field in 
the locale of the bubble were very low, or if 
the local ISM were of low ionization. 

To be yet more quantitative, it is possible to 
relate a measured value of V'max from the bubble 
contribution to the length scale Zb- Equation (Tl8|) 



gives the solution for ^bub in the azimuthal field 
case. Figure E] plots the value of |y max | as a func- 
tion of x. In other words holding x fixed, and 
measuring the polarization PA in a direction per- 
pendicular to the axis of antisymmetry, Figure |S] 
provides the location (x,y max ). 

The scaling for the PA rotation increases 
rapidly with decreasing value of the impact pa- 
rameter, w. Consequently, the strongest measur- 
able effect of the bubble for Faraday rotation will 
be around the center of the projected bubble. At 
these locations one has 9 ~ 0, for which y max is 
easily derived to be 

2/max ~ X/V2, (28) 

which is indicated by the blue line in Figure|5l The 
value of V max is also straightforwardly derivable as 

/lr 2 sini G50 • (29) 

where ^maJl is measured from the data corre- 
sponding to a selection of x on the PA map. From 
these measures one can determine Z\>j sini from 
the data. In the absence of any other information, 
there is an inclination ambiguity. Still, an upper 
limit to Zb is obtained from assuming an edge-on 
view of the system, in which case a lower limit to 
the stellar surface magnetic field strength is also 
obtained. 

The overall scale of the PA rotation through a 
stellar bubble is set by a combination of stellar 
and wind parameters and the size of the bubble 
itself (which, of course, is related to its age). For 
the azimuthal field case, which is the one most 
likely to be of observational significance, we have 
already that B v = B c s (i?»/r). For the number 
density of electrons in a spherical stellar wind, we 
have been using n = n w (i?»/r) 2 . The density 
scale constant is given by 



2 The reader might be concerned that the value of i/j 

max ap- 
pears to diverge as x — > 0, for which y max — > as well. 
First, the apparent divergence does not take place owing 
to the finite size of the star. Second, although large values 
of V'max can be achieved, the polarization PA suffers a 180° 
degeneracy. Many rotations in i/> over a small observed an- 
gular scale rapidly leads to polarimetric cancellation. For 
expected star, wind, and bubble parameters, the vast ma- 
jority of the bubble will like have at best a modest value of 
ip, except perhaps at quite long wavelengths. 
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M 



47r n e m H i? 2 ' 



(30) 



where M is the mass-loss rate, fi c is the mean 
molecular weight per free electron, and Woo is the 
wind terminal speed. Then the scale constant for 
Faraday rotation becomes 



6.0° x 



_Bcf£_\ ( R. ) ( 
. 10 GJ \10Q R e J \ 



10" 6 M e yr- 



- 100 km s 



i) (l^c) 



A 



30 cm 



(31) 



This scale parameter is found to be larger for stars 
with stronger surface fields, larger size, and higher 
mass-loss winds; it is smaller for faster winds and 
larger (older) bubbles. It also has the standard 
quadratic dependence on wavelength. 

Using the preceding expressions, one can derive 
a ratio of the stellar radius to the Faraday rotation 
length scale for the bubble, Zb, to be 



R* 

Zb 



n (i0°g) ( 

4 x 10 11 x ^ 



10- 6 M yr- 



. 100 km s 



"T") (lOO fl ) 



A 



30 cm 



(32) 



The scale of this ratio is huge and thus warrants 
a comment of interpretation. Implicit is that the 
length scale for Faraday rotation zt> has been eval- 
uated for the magnetic field and wind density near 
the stellar surface, which are enormous compared 
to interstellar conditions. As a result, Zb is driven 
to incredibly small values on the order of 1 cm 
and less. However, expected values of Faraday ro- 
tation across the bubble are much smaller than at 
the scale of the star. Defining an associated scale 
z' h related to the dimensions of the bubble instead 
of the star gives {R/z' h ) = {R*/z h ) x (R*/R) 2 , 
which when using the parameterization of equa- 
tion (|32|) is of order 10 -2 . In other words the 
length scale associated with rotation of the po- 
larization PA through 180° is roughly 100 times 
larger than the radius of the bubble. 



4. Conclusions 

The goal of this study has been to develop in- 
sight into the possibilities of using Faraday ro- 
tation to probe stellar magnetism in wind-blown 
bubbles. Adopting spherical symmetry and con- 
sidering two field topologies of azimuthal and split 
monopole fields, analytic and semi-analytic results 
were presented in the form of PA maps arising 
from Faraday rotation across a stellar wind bub- 
ble. The field strength for a split monopole de- 
clines too fast to have detectable effects for Fara- 
day rotation. Instead, only the azimuthal field 
scenario is likely to lead to detectable signals, 
for which the key result derived here is the an- 
tisymmetric morphology for the PA rotations (or 
equivalently, the RMs) across the bubble in addi- 
tion to a simple sinz scaling of the amplitude of 
PA rotation. Inclusion of a central spherical re- 
gion of zero magnetic field leads to prominent and 
antisymmetric PA rotations in arc-like structures 
like those observed in SNR G296. 5+10.0 (Harvey- 
Smith et al. 2010). 

To date, studies of Faraday rotation effects 
across Hn regions appear consistent with vari- 
ations in the interstellar magnetic field (e.g., 
Harvey-Smith, Madsen, & Gaensler 2011). Clearly 
there may be only a restricted subset of bubbles 
for which our models will have applications. More- 
over, dispersion in the interstellar RM may present 
challenges in detection and interpretation of the 
effects we have described even when present, as 
for example variations in the interstellar RM can 
have typical amplitudes of order 10 rad/m 2 (Mao 
et al. 2010). We are considering strategies for in- 
terpreting the interaction of interstellar magnetic 
fields with stellar bubbles using Faraday rota- 
tion effects. Given recent results in relation to 
SNR G296.5+10.0, the PN cases (Ransom et al. 
2008, 2010), and the Rosette Nebula (Savage et al. 
2012), Faraday rotation offers great promise as a 
means of using stellar bubbles to discern the prop- 
erties of stellar and/or interstellar magnetic fields. 
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lems involving Faraday rotation and also to an 
anonymous referee for comments that have im- 
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A. Appendix: Generalized Solution for a Stellar Wind with an Azimuthal Magnetic Field 



For an azimuthal magnetic field of the form B v oc 1/r, there is a general solution to the PA distribution 
arising from Faraday rotation for a spherically symmetric envelope with an electron density that is a power- 
law in radius, n G oc r _m , with m the power-law exponent. Consider a bubble of outer radius R2 and inner 
radius R\. For r < R±, the interior region makes no contribution to the Faraday rotation. Using 8 as the 
variable of integration, with tan.0 = w/z, a sightline enters the outer bubble edge at #2 and enters the inner 
region at 9±. Then the solution for the PA rotation arising from the bubble contribution is 



Although the preceding integral is mathematically valid, it may not be physically plausible for arbitrary 
values of m. The selection of B v oc r 1 essentially assumes a frozen- in field that is dragged out with the 
wind plasma in a constant expansion flow (e.g., Weber & Davis 1967), implying an inverse square-law density. 
However, m^2 could represent changes in the ionization state of the gas with distance from the star. 

The integration of eqn (|A1[> takes into account the front-back symmetry of the situation by integrating 
only over the front hemisphere facing the observer. This gives rise to the factor of 2 appearing in the first 
fraction. The integration limits also allow for an interior central cavity of zero field (or, alternatively, zero 
ionization) to the wind-blown bubble as was discussed in the swept-up wind case of section |2"U1 

With or without a central cavity, the upper limit to the integral is always given by 



for a bubble of outer radius i?2- The presence of a cavity only affects the lower limit to the integral. If there 
is no cavity, then 9% = ir/2 for all sightlines. However, within an interior cavity of radius Ri, the lower limit 
of 9\ = 7r/2 holds only for sightlines that fail to intercept the cavity region, with w > R\. Sightlines with 
w < R\ pass through the cavity; the upper limit to the integration now becomes 




(Al) 



sin 82 = Z&/R2, 



(A2) 



sin#i = vj / R\. 



(A3) 
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Fig. 1. — Illustration of the geometry for evaluating PA changes across spherical bubbles because of Faraday 
rotation. The observer is located off to the right. Linearly polarized radiation moving left to right passes 
through the bubble. Along an arbitrary sightline of impact parameter w, the path length through the bubble 
will be given by z — (—zo) = 2z = 2y 'R 2 — ru 2 , where R is the radius of the bubble. Upon emerging, the 
orientation of the linear polarization against the sky will have rotated in a manner that depends on the 
magnetic field and electron density along the path. 
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Fig. 2. — A false color image with contours overlaid indicating the amplitude of the PA changes that would 
be observed across a wind-blown bubble with an azimuthal magnetic field. The plot is for the bubble 
contribution only. The green and blue colors are negative PA rotations; yellow and red are for positive 
values. The plot space is for the plane of the sky as y versus x, with the coordinates normalized to the radius 
of the bubble R. Note the antisymmetric morphology, and that PA — 0° along the y = axis. Further, this 
morphology is valid for all viewing inclinations, with the amplitude of the PA rotation scaling linearly with 
sini. 
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Fig. 3. — Similar to Fig. [2 this PA map is for a split monopole field. In contrast to the azimuthal field 
case, the PA rotations are all of the same sense, either clockwise or counterclockwise (depending on the 
polarity of the field axis directed toward the observer). Sightlines that fail to intercept the magnetic equator 
when \x\ > Rcosi (see text) make no contribution to a net Faraday rotation, leading to the appearance of 
"truncated" hemispheres; the example shown is for i = 60°. 
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Fig. 5. — (a) Left is a PA map for a bubble with an azimuthal field, like Fig [3j but now with interstellar 
polarization effect included (see text), (b) Right shows how the polarization PA changes as a function of 
y along a horizontal line with x = 0.1R. Black is the full profile. Green signifies the term StpiSM- The 
horizontal blue line is the half-point between the two local maxima. Then the vertical blue arrows represent 
the full change in PA from one side of the bubble to the other because of ipbub(x,y). 
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Fig. 6. — Solution for the location of y m a.*(x) where ^ ma x is achieved in the case of a wind-blown bubble with 
an azimuthal magnetic field (but no cavity). The black line is the numerical solution. For x/R <C 1, y m ax/-R 
is also small, indicating that ipmax occurs near the projected center of the bubble. In this limit y max ~ x/V%, 
as shown by the blue line, which has the correct slope but is shifted upwards slightly for ease of viewing. 
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Fig. 7. — Similar to Figure |6l but now for the case of a swept-up wind. In this example the field-free cavity 
has a radius R\ — 0.75i?2- Solid shows the location of y max (x). Also plotted is the corresponding value 
ro mH = \J x 2 + y2i ax in red. The latter is seen to be flat for \x\ < R\ indicating that the location of the peak 
value V'max follows a circular arc for this range of a:- values. 
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